Chemical control of polar behavior in bicomponent short-period superlattices 
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Using first-principles density functional calculations, we study the interplay of ferroelectricity and 
polar discontinuities in a range of 1-1 oxide superlattices, built out of ferroelectric and paraelectric 
components. Studies have been carried out for a varied choice of chemical composition of the 
components. We find that, when polar interfaces are present, the polar discontinuities induce off- 
centric movements in the ferroelectric layers, even though the ferroelectric is only one unit cell 
thick. The distortions yield non-switchable polarizations, with magnitudes comparable to those of 
the corresponding bulk ferroelectrics. In contrast, in superlattices with no polar discontinuity at the 
interfaces, the off-centric movements in the ferroelectric layer are usually suppressed. The details 
of the behavior and functional properties are, however, found to be sensitive to epitaxial strain, 
rotational instabilities and second-order Jahn- Teller activity, and are therefore strongly influenced 
by the chemical composition of the paraelectric layer. 

PACS numbers: 73.20.-r, 77.84.-s, 71.15.Nc 



I. INTRODUCTION 



Superlattices formed by layer-by-layer epitaxial growth 
of perovskite-based oxide materials are currently a sub- 
ject of intense research, because of their promising tech- 
nological applications as well as fundamental scientific 
interest^. In AB0 3 perovskites, the A +2 B +4 3 (II- 
IV) structures consist of (100) layers of formally charge 
neutral AO and B0 2 , while A+ 3 B+ 3 3 (III-III) or 
A +1 B +5 3 (I-V) structures have charged planes, com- 
posed of +1 AO and -1 BO2 layers or -1 AO and +1 
BO2 layers, respectively. By stacking two perovskite 
layers from different charge families along the [001] di- 
rection, one obtains a polar discontinuity at the inter- 
face. Such polar discontinuities have been reported to 
lead to nontrivial local structural and electronic ground 
states^^, which are often not present in the parent bulk 
compoundgSGnil. Investigating the properties of these 
"exotic" local phases has been an increasingly active area 
of research in the past few years, particularly following 
a 2004 reportP of a conducting quasi-two dimensional 
electron gas (2DEG) at the interface between two wide- 
band insulators, LaA10 3 (LAO) and SrTi0 3 (STO). The 
measured mobility and carrier density of the LAO/STO 
interface are an order of magnitude larger than those 
in analogous semiconductor-based systems^. Further- 
more, magnetism,^ superconductivity^ and a rich elec- 
tronic phase diagram^ have also been reported for this 
same system. These fascinating and unexpected effects 
have generated strong excitement, and an intense effort 
is currently devoted to better understanding the fun- 
damental mechanisms of charge compensation at polar 
oxide-oxide interfaces. Parallel to this thrust, from the 
materials-design point of view, it is also important to 
identify new compounds, artificial superlattices or inter- 
faces that might display similar (or possibly even more 
striking) behavior. 



A system can respond in several different ways to avoid 
a so-called polar catastroph^- a divergence in the poten- 
tial caused by the polar discontinuity - at such an inter- 
face between two charge-mismatched perovskites. Com- 
pensation by free carriers was proposed in Ref. [31 and 
is consistent with the observed conductivity at the in- 
terfaces. Other likely possibilities are direct ionic charge 
compensation through mixed valency of the B cation, ion 
intermixing or oxygen vacancies at the interface^, or po- 
lar distortions at the interface. These can be "induced" 
by the polar discontinuity if both materials are paraelec- 
tric (PE) in their bulk phase^, or "natural" if one or both 
components in the superlattice is ferroelectric (FE)p2l. 

The mechanisms underlying induced and natural polar 
distortions can be readily understood in terms of classical 
electrostatics and the modern theory of polarization^, 
whenever the relevant layers in the superlattice are at 
least three or four unit cells thick. In particular, the 
charge mismatch can be interpreted as a polarization dis- 
continuity, which produces macroscopic electric fields in 
one or both components, because the normal component 
of the electric displacement field, D = £ + 4ttP, is pre- 
served at a coherent insulating interface. For smaller 
thicknesses, macroscopic concepts lose their meaning, as 
each film is too thin to identify a well-defined local value 
of the electric field or the polarization. Thus, it remains 
an important question whether the above concepts are 
still valid when the layers in a superlattice are made thin- 
ner and thinner, down to the ultra-thin extreme limit of 
a single unit cell, or whether new phenomena arise that 
radically alter the physics. It is the main scope of this 
work to answer this question by studying some computer 
designed examples of 1-1 superlattices, constructed out of 
ferroelectric and paraelectric layers, with different com- 
position and choices of polar/non-polar combinations. 
We stress that this question has some points of contact 
with, but is largely unrelated to the better-known prob- 
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lem of the critical thickness for ferroelectricity^ since the 
present problem deals with non-switchable polar distor- 
tions driven by a polar discontinuity rather than switch- 
able ferroelectric behavior. 

In addition to issues related to the presence or absence 
of a polar discontinuity, other factors can influence or 
complicate the properties of the 1-1 bicomponent super- 
lattices. In particular, structural instabilities such as the 
rotation and tilting of oxygen octahedra, the presence of 
misfit strain, and the presence of d° ions - associated with 
a tendency to ferroelectricity - in otherwise paraelectric 
layers, are likely to strongly affect the behavior. It is 
therefore of interest to study how these additional factors 
may alter the properties of the short-period superlattices 
compared to those expected solely out of consideration 
of the presence or absence of a polar discontinuity. 

To separate the role of electrostatics from the role 
of the specific choice of cations within a given charge 
family, we investigate bicomponent superlattices with al- 
ternating III-III/II-IV and I-V/II-IV perovskite layers 
(with interfacial polar discontinuities), as well as oth- 
ers with alternating II-IV/II-IV and I-V/I-V perovskite 
layers (without polar discontinuities), both for a variety 
of chemical constituents. We find that our 1-1 structures 
behave in a way which is qualitatively analogous to the 
longer-period superlattices previously investigated in the 
literature. In particular, the systems containing polar 
discontinuities have strongly broken inversion symmetry 
and large polarizations, but are not switchable. The mag- 
nitudes of the polar distortions are as large as those in 
the corresponding bulk FE materials. The superlattices 
formed out of PE and FE layers with no potential discon- 
tinuity at the interface, on the other hand, usually have 
their off-centric displacements suppressed. We do, how- 
ever, find examples of superlattices without polar mis- 
match that show marked ferroelectric instability, driven 
by either suppression of rotational instability or due to 
the presence of d° ions in the PE layer. 

The remainder of the paper is organized as follows. 
Section II includes the methodology and the details of the 
computations carried out in this study. In Section III we 
present our results and discussion in several sub-sections 
dealing respectively with the bulk properties of the con- 
stituent materials, the calculated structural properties 
of the superlattices (with and without polar discontinu- 
ities), and the calculated functional properties including 
static dielectric constants. Section IV provides the sum- 
mary. 



II. METHODOLOGY AND COMPUTATIONAL 
DETAILS 

Our calculations are performed using the plane wave 
based pseudo poten tial framework of density functional 
theory (DFTjpSECD as implemen ted in the Vienna ab ini- 
tio Simulation Package (VASPJpEl We choose the lo- 
cal density approximation (LDA) rather than generalized 



TABLE I: Representative bicomponent superlattices com- 
posed of paraelectric (PE) and ferroelectric (FE) components. 





PE 


FE 


Polar Interface? 


Mismatch (%) 


LAO/PTO 


LaA10 3 


PbTiOa 


YES 


0.99 


KSO/PTO 


KSbOa 


PbTiOs 


YES 


1.94 


KSO/KNO 


KSbOs 


KNb0 3 


NO 


0.10 


SSO/PTO 


SrSn0 3 


PbTiOs 


NO 


4.40 



gradient approximation (GGA) to describe the exchange- 
correlation functional due to the later's systematic over- 
estimation of polar behavior^. The optimized geome- 
tries are obtained by full relaxation of the atomic posi- 
tions and out-of-plane lattice constants (c), with the in- 
plane lattice constants (a) fixed to the LDA bulk lattice 
constant of the ferroelectric component in the tetrago- 
nal phase. The positions of the ions are relaxed toward 
equilibrium until the Hellmann-Feynman forces become 
less than 0.001 eV/A. A 4x4x4 Monkhorst-Pack k-point 
mesrPil and 450 eV plane-wave cut-off give good conver- 
gence of the computed ground state properties. 

The Born effective charges, phonon frequencies and 
dielectric tensors are obtained from linear response 
calculations^!, using a variational formalism^ of the den- 
sity functional perturbation theory, as implemented in 
the ABINIT code. We calculate the electronic contri- 
bution to the polarization as a Berry phase using the 
method developed by King-Smith and VanderbiltpS, and 
extract the ionic contribution by summing the product 
of the position of each ion with its pseudo-charge. 

The first part of our study concerns superlattices 
formed by alternating FE and PE single simple per- 
ovskite unit cells in the z direction, and constraining 
them to be tetragonal in the x — y plane. In a sec- 
ond stage, we consider structures formed by alternating 
2x2x2 PE and FE cubic perovskite unit cells to allow 
rotation and tilting of the BC>6 octahedra. 



III. RESULTS AND DISCUSSION 

We consider the four representative systems listed in 
Table |TJ All combinations contain one material which 
is ferroelectric in the bulk, and one which is paraelec- 
tric. Two of the systems (LAO/PTO and KSO/PTO) 
have polar discontinuities at the interface whereas the 
other two (KSO/KNO and SSO/PTO) do not. We 
deliberately choose paraelectric components without d° 
cations, because d° cations have a strong tendency to 
off-center through the so-called second-order Jahn- Teller 
effect. This often results in ferroelectricity, and in mate- 
rials with paraelectric ground states it can cause incipient 
or quantum ferroelectricity^, which we choose to avoid 
here. We begin by briefly summarizing the relevant bulk 
properties of the various constituents. 
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FIG. 1: (Color online) Upper panel: Polar catastrophe for atomically abrupt (001) interfaces between LAO and PTO, KTO 
and PTO, KSO and KNO, SSO and PTO (left to right). The charges of each AO and B0 2 layer are taken from the formal 
charges of the constituent ions, after Ref. |1 Middle panel: Optimized structures for LAO/PTO, KTO/PTO, KSO/KNO 
and SSO/PTO, projected onto the xz plane. The large spheres indicate the cations at various layers, while the small spheres 
indicate oxygens. Inequivalent oxygens are indicated. Lower panel: Cation displacements along the z axis, with respect to the 
centers of their oxygen cages. 



1. Bulk properties 

Both PTO and KNO are well-known ferroelectric ma- 
terials with the cubic perovskite structure at high tem- 
perature. PTO undergoes a phase transition to a tetrag- 
onal phase below 766 K, with the polarization along the 
[001] direction. In the tetragonal phase the experimental 
in-plane (a) and out-of plane (c) lattice parameters are 
3.90 A and 4.14 A respectively^ and the ferroelectric po- 
larization is 59 ^C/cm -2 ; our calculated corresponding 
LDA values are 3.85 A and 4.05 A, and 80.45 ^iC/cm" 2 
respectively. Bulk KNO is a rhombohedral ferroelectric 
with a polarization of 42 ^tC/cm~ 2 along [111] below its 
Curie temperature^ of 230 K. Earlier theoretical work 
reported stabilization of the tetragonal phase with com- 
pressive strairP^l for tetragonal KNO we obtain in- and 
out-of-plane lattice constants of 3.92 A and 4.00 A and 
a polarization of 26 fiC/cm 2 along [001]. 

LAO is a wide band gap insulator with strong al- 
ternating rotations of the oxygen octahedra around the 
[111] direction leading to rhombohedral symmetry^. Re- 
cent density functional calculations showed that epitax- 
ial strain causes a polar instability which competes with 
these non-polar oxygen rotations^. The lattice parame- 
ter for the cubic structure without oxygen rotations is 
3.81 A^. KSO is not known experimentally. Com- 
putationally we obtain a cubic LDA optimized lattice 
constant of 3.92 A. At room temperature SSO crystal- 
lizes in the centrosymmetric orthorhombic space group 



Pbnm&l, with lattice parameters a — 5.71 A, b = 5.71 
A and c = 8.06 A. The structure is characterized by a 
classic GdFe03 tilt pattern of the SnOe octahedra. A 
structural transition from the ground-state orthorhom- 
bic structure to a high temperature tetragonal 1 4/ mem 
phase has been founcP^ driven by an order-disorder oc- 
tahedral tilting transition. 



2. Structural properties 

We begin by optimizing the geometries of our four rep- 
resentative superlattices. We construct the superlattices 
by alternating FE and PE layers along the [001] direc- 
tion. We impose an in-plane periodicity of one simple 
cubic perovskite unit cell: This choice imposes an overall 
tetragonal symmetry and prohibits rotations and tiltings 
of the oxygen octahedra. We investigate the effects of re- 
laxing this constraint later. By fixing the in-plane lattice 
constant to the LDA value for the FE layer we induce 
a small tensile strain in the PE layer for LAO on PTO, 
small compressive strains for KSO on PTO or KNO, and 
a large compressive strain for SSO on PTO (values in Ta- 
ble [I]). We then fully relax the out-of-plane lattice param- 
eter and the ionic positions within the tetragonal sym- 
metry. All superlattices are found to be insulating within 
the LDA with band gaps of ^2 eV, allowing us to directly 
calculate their ferroelectric polarization using the Berry 
phase formalisnPll. 
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Both superlattices containing polar discontinuities are 
strongly polar, with polarizations of -50.87 /xC/cm 2 for 
LAO/PTO and +57.74 /iC/cm 2 for KSO/PTO; the op- 
posite signs in the two cases reflect the opposite ori- 
entation of the polar discontinuity (III-III /II-IV in the 
first case and I-V/II-IV in the second). Neither case 
is ferroelectric, because the polarization orientations are 
fixed by the polar discontinuity and the polarizations are 
not switchable^. The calculated structures are shown 
in Fig.l. In LAO/PTO, the Pb and Ti ions displace 
from the centers of their coordination polyhedra by w0.25 
A and 0.15 A respectively along the positive crystallo- 
graphic z axis, that is towards the PbO/A102 interface. 
In KSO/PTO, Pb and Ti off-center by 0.29 A and 0.18 
A respectively along the negative z axis, away from the 
PbO/Sb02 interface. The magnitudes of these displace- 
ments are similar to those in bulk PTO, which are 0.41 
A and 0.26 A for Pb and Ti respectively. Note that the 
off-centric movements in the LaO/KO and A102/Sb02 
layers are tiny, with the movements of K and Sb being in 
the same direction as Pb and Ti, and La and Al moving in 
the opposite direction to each other. We, however, do not 
attach much significance to these small displacements. 

The situation for the superlattices without polar dis- 
continuities is more complicated. For KSO/KNO we in- 
deed obtain a negligible polarization of 0.02 /LtC/cm 2 , re- 
sulting from the tiny displacements of K and Nb from 
the centers of their coordination polyhedrEp^. Note that 
the KSO/KNO superlattice in its paraelectric phase has 
an inversion center. However we do not impose the in- 
version symmetry during the relaxation procedure, since 
our aim is to find out the existence or non-existence of 
the polar behavior. The tiny polarization obtained for 
relaxed geometry of KSO/KNO is within the accuracy 
limit of our calculations. This leads us to conclude that 
our obtained geometry is in fact centrosymmetric, as has 
been checked in terms of the energy difference between 
the fully centrosymmetric case and the relaxed struc- 
ture, which turned out to be less than 0.01 meV. For 
SSO/PTO, however, we obtain a rather large polariza- 
tion of ~ -66 ^C/cm 2 , with a cooperative displacement 
of all cations (Sr, Sn, Pb, Ti). Interestingly, the polar- 
ization is found to be switchable: In Fig. 2 we show the 
calculated total energy, £ as a function of polarization 
P obtained by calculating polarizations and total ener- 
gies of geometries with fixed values of off-centric move- 
ments. The fixed values of off-centric movements are 
obtained by interpolating between off-centric movements 
of ions corresponding to two minima of the E versus P 
curve. The E(P) curve shows the characteristic double 
well structure typical of switchable ferroelectrics, except 
that in this case it is strongly asymmetric as a result 
of the additional inversion-center lifting that is built in 
to the structure through the layering. The switchable 
behavior indicates that the polarization arises from an 
origin other than the polar discontinuity alone, as the 
polar discontinuity can only drive the polarization in one 
direction. Note that it is also likely distinct from the be- 
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FIG. 2: Potential energy curve of the total polarization for 
SSO/PTO superstructure for lxlxl/lxlxl geometry (see 
text for details). 



havior in BaTiOa/SrTiOs superlattices^ where the PE 
component (STO) is a quantum ferroelectric. In the next 
section we explore the origin of this behavior. 



3. Origin of polarization at non-polar interfaces 

Our first clue as to the origin of the unexpected po- 
larization in the SSO/PTO superlattice lies in the large 
(4.4%) mismatch in lattice constants between the two 
constituents in this system. To test the role of the misfit 
strain in driving polarization, we performed an analogous 
set of calculations for SrSn0 3 /BaTi0 3 (SSO/BTO), in 
which the misfit strain ( estima ted using the experimen- 
tal bulk lattice constants^ 1221 ) is only 0.78 %. Indeed, 
we find that the calculated polarization for SSO/BTO is 
considerably reduced (to -7.76 /iC/cm 2 ) compared with 
SSO/PTO, and the corresponding off-center displace- 
ments are 0.02 A, 0.06 A, 0.06 A and 0.02 A for Sr, 
Sn, Ba and Ti respectively. We have further carried out 
calculations for bulk SrSn03 in tetragonal geometry, dis- 
abling the rotation and tilt of the SnOg octahedra, under 
4.4 % and 0.78 % compressive strains (strains identical to 
SSO/PTO and SSO/BTO systems respectively). These 
calculations resulted in large off-centric displacements of 
Sr 2+ and Sn 4+ ions for 4.4 % compression strain and neg- 
ligible displacements for 0.78 % strain value, illustrating 
once more the essential role of the strain in determining 
the polarization. 

While it is an already established resullP^ that strain 
plays a crucial role in determining polarization and 
off-centering in superlattices, there has also been re- 
cent discussion of competition between tilts and rota- 
tions of the oxygen octahedra and ferroelectric displace- 
ments. For example it has been shown computationally^ 
that LaA103 shows ferroelectric behavior under epitax- 
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ial strain when the octahedral rotations are suppressed, 
while ferroelectricity is suppressed when those rotations 
are allowed to relajp^. Similar physics led to predic- 
tions of ferroelectricity in A-site alloy perovskites^, in 
which the polar behavior is obtained by destabilizing the 
cooperative rotational instability through A-site disor- 
der. And in Ref. I39|, polar behavior was engineered the- 
oretically in CaMn0 3 by suppressing the orthorhombic 
tilting instability through strain and chemical substitu- 
tion. The possibility of competing rotational instabil- 
ity governed by zone-boundary phonons with ferroelec- 
tric modes should therefore be carefully analyzed before 
drawing a definitive conclusion on a specific system. 

To check for this possibility in our superlattices, we in- 
creased the size of our supercell to 2 x 2 x 2 simple per- 
ovskite unit cubes in order to allow for oxygen rotations 
and tiltings and re-compute the lowest energy structures. 
Indeed we find that, while for LAO/PTO, KSO/PTO 
and KSO/KNO our conclusions remain unchanged, for 
SSO/PTO the system adopts a centrosymmetric struc- 
ture with negligible off-centering when rotational distor- 
tions are allowed. Bulk SrSn0 3 calculations allowing for 
the octahedral rotation and tilt, as expected, also yielded 
a centrosymmetric structure. We find applied compres- 
sive strain increases the rotation and decreases the tilt 
in agreement with the findings in Rcf 29 > This reinforces 
the earlier suggestions in the literature, and implies that 
if octahedral rotation at the SSO layer can be disabled, 
which might be possible in ultra-thin films using strain or 
doping, then SSO/PTO superlattices could have switch- 
able polar properties. 



4- Ferroelectric tendency in the paraelectric layer 

In all of the model superlattices studied above, the 
paraelectric material was chosen to be one that shows no 
tendency to ferroelectric off-centering in its bulk phase. 
In this section we study the effect of substituting it with 
an "incipient ferroelectric" or "quantum ferroelectric" , 
i.e. a material that is still PE in the bulk, but very close 
to a ferroelectric phase transition. 

Specifically we compare our results reported 
above for KSO/KNO with analogous calculations 
for KTa0 3 /KNb0 3 1-1 superlattices. We choose KTa0 3 
to compare with KSb0 3 because the former has a 
so-called Second Order Jahn- Teller (SOJT) active^ ion 
with a d° electron configuration which is associated with 
a tendency to off-center. Indeed bulk KTO is a quantum 
ferroelectric with ferroelectricity suppressed by quantum 
fluctuations^. Ultrashort-period (KTO)™/KNO super- 
lattices have been previously investigated^ in terms of 
ab-initio studies. As in Ref. f4"Tl for m = 1, we find that 
KTO/KNO is strongly ferroelectric, with a switchable 
polarization of 16.46 /xC/cm 2 and atomic displacement 
at the Ta site of 0.14 A. 

The difference in physical behavior between 
KTO/KNO and KSO/KNO must arise from the 



chemical difference between KTO and KSO. To further 
understand this, we next compute the Born effective 
charges, Z* for KTO and KSO, using the c/a ratios 
that we obtained for the superlattice geometries. The 
Ta Z* in KTO is highly anomalous (8.78 instead of the 
nominal valence of 5), consistent with the proximity 
of the d° KTO system to ferroelectricity, while the Sb 
Z* in KSO is 5.4, close to the nominal value of 5. The 
strong ferroelectric instability in KNO is therefore able 
to induce a ferroelectric instability in the proximally 
ferroelectric KTO layer; KSO is less polarizable and 
therefore resists polar distortion. 

We also compute the properties of the KTa0 3 /PbTi0 3 
superlattice in order to study the behavior of a ferroelec- 
trically active paraelectric in a superlattice with a polar 
discontinuity. Here KTO is under 3.51 % compressive 
strain. We find a non-switchable polarization of 69.72 
fiC/cm 2 and displacements of 0.39 A, 0.26 A, 0.06 A and 
0.10 A for Pb, Ti, K and Ta respectively; all of these 
quantities are larger than the corresponding values for 
LAO/PTO and KSO/PTO, expected because in this case 
SOJT activity is additive with polar discontinuity behav- 
ior. 



5. Static Dielectric constant 

Finally, for completeness, we compute the dielectric 
properties of the polar PE-FE superlattices driven by 
polarization discontinuity, specifically LAO/PTO and 
KSO/PTO. Our motivation lies in the fact that the di- 
electric constants typically hold signatures of polar in- 
stability of soft modes^ and are also of technological 
importance. 

The a/3 component of the static dielectric constant is 
given by, 



£uf3 — { e<X )al3 + 



-in 



E 



(i) 



where e a is the electronic contribution to the dielectric 
tensor and the second term is the sum of the contribu- 
tions from each polar phonon. /i labels the phonon modes 
with frequency u; M , Z* Q is the mode effective charge cor- 
responding to mode /i in the Cartesian direction a and 
51o is the volume of the primitive unit cell. The mode 
effective charge in the a direction for a given mode n is 
related to the eigendisplacement involving ion k 

and the Born effective charge tensor, Z* by, 



1/2 • 



(2) 



Electronic dielectric constant: The results for the elec- 
tronic dielectric tensor are shown in Table-II. The elec- 
tronic contributions are of almost the same magnitude 
for LAO/PTO and KSO/PTO superlattices as expected 
from the similarity of their LDA band gaps. 
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TABLE II: Electronic and lattice contribution of the static di- 
electric constant for LAO/PTO and KSO/PTO superlattices. 





Electronic contribution 
LAO/PTO KSO/PTO 


Lattice contribution 
LAO/PTO KSO/PTO 


€xx 

e yy 


6.39 5.93 
6.39 5.93 
6.14 5.34 


137.19 90.12 
137.19 90.12 
43.09 24.99 



TABLE III: Phonon frequencies in cm -1 , corresponding mode 
effective charges in unit of je| and the contribution to static 
dielectric constant. 



LAO/PTO 


KSO/PTO 


oj 


7* e zz 

^ fJ,,ZZ 


LU 


7* 

^ fl,ZZ 


zz 


115.19 


0.82 21.95 


125.83 


0.57 


8.32 


181.17 


0.70 6.50 


194.37 


0.60 


4.03 


302.57 


1.24 7.33 


209.94 


1.06 


10.06 


363.04 


0.74 1.83 


284.69 


0.36 


0.68 


564.96 


0.60 0.49 


458.66 


0.43 


0.37 


681.35 


1.01 0.95 


639.09 


1.08 


1.20 


757.08 


0.27 0.05 


801.07 


0.84 


0.46 



Lattice dielectric constant: Our calculated lattice con- 
tributions to the dielectric tensor are listed in Table-II. 
We find that the lattice contribution to the static di- 
electric constant is quite different between LAO/PTO 
and KSO/PTO superlattices. To understand the origin 
of the difference, we next carry out detailed analysis of 
the phonon frequencies and mode effective charges for 
the phonon modes which contribute to e zz ; that is the 
relevant quantity for superlattices grown along the [001] 
direction with significant off-centric movements along the 
z-axis. A large lattice dielectric response can be the re- 
sult of the presence of one or more very low frequency 
polar phonons and/or anomalously large mode effective 
charges. Our analysis, as listed in Table-Ill, shows that 
for the LAO/PTO superlattice, the lattice dielectric con- 
stant originates essentially from the dominant contribu- 
tion of the lowest mode at 115 cm -1 ; this in turn arises 
from the combined effect of large mode effective charge 
and low frequency of the phonon mode. The contribu- 
tions of different ions to the mode effective charge are 
found to be of comparable magnitudes. 



In the case of the KSO/PTO superlattice, on the other 
hand, the dominant contributions come from more than 
one phonon mode. The major contributions arise from 
the phonon mode at 209.94 cm -1 , which has a giant 
mode effective charge, and that at 125.83 cm' 1 , which 
has low frequency. The giant mode effective charge for 
the phonon mode at 209.94 cm -1 comes primarily from 
the Ti +4 and Sb +5 ions and the planar oxygens situated 
at the Sb02 layer. Thus, the temperature dependence of 
dielectric response of KSO/PTO superlattice is expected 
to be weaker than that of the LAO /PTO superlattice. 
IV. SUMMARY 

We have carried out a detailed computational study 
of superlattices formed out of alternating layers of one 
unit cell thick FE and PE components, considering sev- 
eral different ABO3 perovskites. We find that the ten- 
dency of the superlattices to exhibit polar properties de- 
pends strongly on the chemistry of the PE and FE com- 
ponents. In addition to interfacial polar discontinuities at 
the interfaces between the PE and FE layers, the epitax- 
ial strain, rotational instability of the PE layer and the 
presence of d° second order Jahn- Teller active ions in the 
PE layer contribute to determining the polar behavior. 
While the polarization arising from a polar discontinuity 
is generally not switchable, other routes to polar behav- 
ior discussed here, do allow polarization switching. The 
strongest polar response was found for superlattices with 
polar discontinuities and formed out of FE and d° PE 
layers. 
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